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Abstract 

In the present paper, we study the Poincare map associated to a periodic perturba- 
tion, both in space and time, of a hnear Hamiltonian system. The dynamical system 
embodies the essential physics of stellar pulsations and provides a global and qual- 
itative explanation of the chaotic oscillations observed in some stars. We show that 
this map is an area preserving one with an oscillating rotation number function. The 
nonmonotonic property of the rotation number function induced by the triplication 
of the elliptic fixed point is superposed on the nonmonotonic character due to the 
oscillating perturbation. This superposition leads to the co-manifestation of generic 
phenomena such as reconnection and meandering, with the nongeneric scenario of 
creation of vortices. The nonmonotonic property due to the triplication bifurcation 
is shown to be different from that exhibited by the cubic Henon map, which can be 
considered as the prototype of area preserving maps which undergo a triplication 
followed by the twistless bifurcation. Our study exploits the reversibility property of 
the initial system, which induces the time-reversal symmetry of the Poincare map. 
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1 Introduction 



During the past twenty years there has been increasing interest in the study of 
the coexisting regular and chaotic motions in Hamiltonian systems with few 
degrees of freedom. A basic problem of study was the transition to chaos in 
two degrees of freedom Hamiltonian systems which are slight perturbations of 
an integrable system, whose Hamiltonian is nondegenerate. In order to ease 
this task, most of the results on the dynamics of such Hamiltonian systems 
were obtained exploring the behaviour of a Poincare map associated to a 
cross section in a fixed energy surface. A Poincare map is an area preserving 
map (APM) defined on an annulus, that is, depending on {r,6), where r is a 
momentum-like coordinate and 6 is an angular variable. Considering the lift 
of such a map, namely computing the angle without the restriction modulo 
27r, one associates to an orbit (r„, 9n) the rotation number which is the limit: 

p = hm (1) 



if it exists. If p is an irrational number, then the orbit densely fills an invariant 
circle (KAM circle), while an orbit of rational rotation number p/q — with 
p, q coprime integers, ^ > — is a ^-periodic orbit. 

Typical questions of mathematical and physical interest include the persistence 
of the KAM circles after perturbation and determination of the threshold at 
which a circle of a given rotation number breaks up. The main results which 
answer these questions, such as the KAM theorem, the Poincare-Birkhoff 
theorem, and the Moser twist theorem, are based on the validity of the twist 
property of the area preserving map transforming (r, 9) ^ (r', 9'), namely: 

B9' 

^v^O,Vr (2) 



The twist condition implies the monotonic change of the rotation number p 
with r. The APM satisfying this condition is called monotonic twist map or 
simply a twist map. During the last twenty years the monotonic twist APMs 
were studied extensively [1]. Consequently, their dynamics has gained almost 
complete understanding. In particular, this statement is a consequence of the 
late understanding that Hamiltonian systems have universal behaviour which 
can be revealed by studying explicit symplectic maps of the plane. 

Recently, due to the increasing number of physical phenomena that are mod- 
eled by nontwist maps, this class of dynamical systems has finally captured 
the attention of the scientific community. The few studies on nontwist maps 
existent in the literature have revealed an unexpected variety of local and 



2 



global bifurcations: dimerized island chains, periodic orbits collision, separa- 
trix reconnection, vortex islands [2,3,4]. 

In the present paper we study the Poincare map associated to a periodic per- 
turbation, both in space and in time, of a linear Hamiltonian system. The 
dynamical system originates from a model of stellar pulsations previously in- 
troduced in [5] for low mass stars {M/Mq < 8, where Mq is the solar mass) 
and extended to intermediate mass stars {M/Mq G [8, 11]) in [6]. The latter 
work included some preliminary results concerning the nontwist character of 
the map. Here, wc show that this map is an APM with an oscillating rotation 
number function. In spite of the system being a perturbation of a linear Hamil- 
tonian, we argue that due to both the acquired nontwist character and the 
periodic character of the perturbation, the map presents the typical features of 
a generic class of nontwist APMs, such as reconnection and meandering, with 
the nongeneric scenario of creation of vortices. The nonmonotonic property 
due to the triplication bifurcation is shown to be different from that exhibited 
by the cubic Henon map, which can be considered as the prototype of APMs, 
and which undergoes a triplication followed by the twistless bifurcation. Our 
study exploits the reversibility property of the initial system, which induces 
the time-reversal symmetry of the Poincare map. The paper is organized as 
follows. In §2 we introduce our model of stellar pulsations. Before getting into 
the detailed analysis of it, we consider worth recalling the most important 
features of nontwist maps. This is done in §3. In this section we will also 
emphasize the consequences of both the triplication bifurcation and of the 
symmetry properties. Later, in §4 we study the dynamics of our simple model 
of stellar pulsations. Finally in section 5 we summarize our major findings and 
we draw our conclusions. 



2 Description of the model 

Pulsating stars are perfect tests of the theory of stellar evolution since the 
comparison between the theoretical pulsational characteristics (periods, am- 
plitudes, and growth time scales) and the observed ones can help substantially 
in the fine tuning of the stellar models. Low- and intermediate-mass stars are 
prone to experience recurrent thermal instabilities and substantial mass loss 
during certain intervals of their evolution. Moreover, a substantial amount of 
pulsating stars in the Galaxy shows irregular behaviour and this makes them 
interesting targets in the field of dynamics systems. The present model em- 
bodies the essential of stellar oscillations by assuming an oscillatory driving 
originating in the stellar interior in the form of sinusoidal pressure waves. This 
model belongs to the family of the so-called one-zone models which treat the 
star as a rigid core surrounded by a dynamic, homogeneous gas shell. The 
pressure waves propagate through a transition zone characterized by a certain 



3 



transmission coefficient until they hit the stellar envelope and dissipate. No 
back reaction of the outer layers on the inner ones is considered. The dissipa- 
tion of the pressure waves induces a fluctuation in the radius and velocity of 
the outermost layer and this variability makes the object of our study. 

The variation of the interior radius, Rc, around an equilibrium value, Rq, is 
given by Rc = Ro + aRo sin uJcT, where r is the time in years and a and Uc 
are the fractional amplitude and the frequency of the driving, respectively. In 
absence of any driving force, the equation of motion is 

d^R _ AttR^ GM 



where M is the stellar mass and R and P are the radius and pressure in the 
shell of mass m, respectively. For convenience, we introduce the nondimen- 
sional variables x = R/ R^ — 1 and t = cu^t, where R^, is the equilibrium 
stellar radius and 



GM 



Ri 



■■k 



is the characteristic frequency of the star. 

Considering that the additional perturbative acceleration is proportional to 
the driving acceleration with a proportionality coefficient Q we obtain 



(4) 



ret 



and after some algebraic manipulation [6] the final equation of motion reads 
X — y 

(5) 

y — —X — e sin{ujx — ut — auj^^^ sin ut) , 

where the characteristic frequency of the system, uj results from the use of the 
dimensionless time unit and from the assumption that Rq encompasses almost 
the entire stellar mass. It is defined as 

In the action-angle coordinates (r, 9) resulting from the change x — V2rcos9 
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and y — -\/2r sin 9, the Hamiltonian of the system is 

H{r,6) = r — — cos{u!\/2r cos6 — ut — auj^^^ sinoot). (7) 

LO 

The parameters of the system are e, the total driving amphtude, a, the frac- 
tional driving amplitude and uj, the characteristic frequency of the system. 
From the point of view of the parametric study, we take u to be the most 
important parameter fixing the type of star to be studied and its evolution- 
ary status. It is a measure of the ratio between the internal radius where the 
driving is produced and the equilibrium stellar radius. From the astrophysi- 
cal perspective, we are interested in the dynamics of the system from Eq.(5) 
corresponding to c<j ^ 3, more exactly u = 3.0146. This low value of uj is asso- 
ciated with evolved stages of intermediate mass stars, as it results from stellar 
evolutionary models [7,8]. Since we are interested in modeling the evolution 
of this type of stars, we fix u! during our simulations. We have performed a 
thorough parametric study in the space (e,a) as it assures a fine tuning of the 
strength of the perturbation. As in [6], we are interested in e < 100% and 
a < 40%, which allow the perturbative approach. 



3 Nontwist maps 

The prototype of nontwist area preserving maps is the quadratic standard 
map [4,9]: 

r„+i — Tn — k sin 6 

n+l n 

On+1 = On + 27ra; - r^^^ (mod 27r) . 



The twist condition is violated along the curve r = k sin 6. For A; = we get an 
integrable map, whose orbits lie on the circles of constant r. The rotation num- 
ber function is p(r) = a; — r^/(27r). The circle r — has the maximum rotation 
number p^ax = tu. It is called the twistless circle or shear less circle. A slight 
perturbation leads to the persistence of the twistless circle and nearby circles 
having a diophantine rotation number — that is, an irrational number badly 
approximated by rationals [10]. Hence, the nontwist standard map, defined 
above, has an invariant circle of maximum rotation number among the rota- 
tion numbers of the nearby orbits. For a fixed perturbation parameter k and 
for UJ chosen in such a way that the map has a pair of Poincare-Birkhoff chains 
containing p/q-tjpe periodic orbits, as uj varies the two chains approach each 
other and at a threshold uj = ujr{k) the hyperbolic points of the two chains 
enter the twistless circle. This is called the reconnection threshold [12]. At the 
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Fig. 1. The triplication scenario. The phase portrait around the eUiptic point at 
the tripHcation threshold is shown in (b). Observe the change in the position of the 
heteroclinic triangle in (a) and (c). 

reconnection threshold the two distinct hyperbolic orbits, of rotation number 
p/q, are connected by heteroclinic arcs. 

There are many ways in which the monotonic twist condition can be violated. 
The vertical lines 6 = Oq can be mapped into curves having a single extremum 
— quadratic twist [9] — or multiple cxtrema — cubic [11], quartic [12] or 
sinusoidal twist [13]. Moreover, it has been proved [14] that whenever the 
elliptic fixed point of an APM of the plane passes through a triplication, a 
twistless bifurcation occurs or, equivalently, the rotation number as a function 
of the distance from the elliptic fixed point becomes a nonmonotonic function. 
The triplication of the elliptic fixed point of an APM of the plane, occurs at 
the value /i = /xs at which the multipliers of the elliptic point cross the values 
A = exp(±27ri/3). At the triplication threshold fj, = /i^ an unstable period-3 
orbit emerges from the elliptic fixed point in both directions, that is for /i < /i^ 
and for /j, > /i^. As /i increases the elliptic point crosses the 1/3 resonance and 
the triangular shape figure changes side after shrinking to elliptic fixed point 
(Figure 1). 

The symmetry properties arc very useful for the explanation of the behaviour 
exhibited by the system under study. We next recall some notions and results 
concerning time-reversal symmetry of a dynamical system [15]. A system of 
differential equations 

x^F{x),xeW, (9) 

or equivalently, the vector field F is called reversible if there is a smooth 

involution R : R ^ R {R o R = Id) such that DRo F = -F o R, where DR 
is the differential of the map R. R is called the reversor of the vector field 
F. If $t is the flow of the vector fleld F, then the reversibility means that 
^-tR — R^f This is tantamount to say that the reflection of the trajectory 
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Fig. 2. The triplication bifurcation in the cubic Hcnon map. (a) The separatrix 
associated to the hyperboHc period-3 orbit of the map from Eq.(ll), corresponding 
to /L( = —1.13, and some nearby orbits; {h) The same map at the saddle-center 
threshold /Xg = —1.12189. The twistless circle is the circle with three cusps. 



x{i) = ^t^o is also an orbit of the system. A diffeomorphism / : — > is 
called reversible with respect to a smooth involution R of if — RofoR. 
If / is an i?-reversible diffeomorphism then I — f o R is also an involution 
(/ o / = Id), and f = I o R. This factorization of the map / is very useful for 
the study of its dynamical properties. It can be shown that /" is a i?-rcvcrsible 
map, too, V n e Z. Denote by /„ = /" o i?, n e Z, and by r„ the fixed point 
set of the involution /„. r„ are called symmetry lines of the map /. A point in 
the intersection Fj fl Fj. is a periodic point of the map /, whose period divides 
\j — k\. The symmetry lines are transformed by /" into other symmetry 
lines in the following way: 

= r2n+k (10) 



In order to interpret the nontwist-type dynamics which will be encountered 
in our system when studying the triplication of the elliptic fixed point, we 
consider necessary to illustrate the twistless bifurcation in the case of the 
cubic Henon map 

Xn+l = -Vn + + flXn + 0.7 ^^^^ 
Un+l -^n ) 



which has an elliptic fixed point {xq, Uq — Xq) that can be easily found. At 
/i3 = —1.1505 this point undergoes a triplication, while at = —1.12189 
the saddle center collision occurs. In Figure 2, we show that as fi increases 
beyond the elliptic and hyperbolic period-3 orbits collide and disappear in 
a saddle-center bifurcation at /i — /is- 
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Fig. 3. Rotation number p{x, ji) versus x for three values of between the triphcation 

and saddle-center collision thresholds of the Henon cubic map (Eq.ll). The curves 
correspond to /x = —1.135, jj, = —1.13 and ji = jjLs = —1.12189. See text for more 
details. 

In order to better identify the twistless bifurcation — the birth of a twistless 
circle — in Figure 3 we show the rotation number p(a;, /x) as a function of the 
points X lying on the symmetry line Fq := {{x,y)\y — x} of the cubic Henon 
map, which is a reversible map with respect to the reversor R{x,y) = {y,x). 
More precisely, the rotation number is computed at the points x G Fq and 
satisfying x < Xq. Figure 3 shows that between the triplication threshold and 
saddle-center collision threshold of the two period-3 orbits, the twistless circle 
is the circle of rotation number pmax =1/3 (the maximum points from right to 
left correspond to /i — —1.135, /i = —1.13, and fx — fig — —1.12189). Observe 
that as /i increases up to /is, the twistless circle moves away from the elliptic 
fixed point. 

Much more important (and difficult) is to detect the location of the twistless 
circle after the collision-annihilation of the two period-3 orbits. In order to get 
some insight into the manifestation of the twist propertiy beyond the threshold 
/is, we have computed numerically the rotation numbers for points starting 
near the elliptic fixed point of the cubic Henon map and on the symmetry 
line Fq. The rotation number functions p{x,ii) are plotted for fi = iig, and 
some parameter values obtained varying fi in the direction given by the vector 
V = fis — (Figure 4). Note that the cubic Henon map has a twistless circle 
of maximum rotation number among the rotation numbers of nearby orbits. 
Varying the parameter /i as specified above, pmax decreases as a result of the 
collision-annihilation of the orbits [4], and at some value, /it, the rotation 
number becomes decreasing. At this point the map is no longer nontwist. This 
is the reason why fi = fit was called the threshold of the twistless bifurcation 
in [14] where the above presented scenario was first presented although in the 
reverse order. 
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Fig. 4. Rotation number p{x,fi) of the points {x,y) on the symmetry line Tq of the 
map where is the Henon map. From top to bottom: /i = —1.12189, fi = —1.12, 
H = —1.115, = —1.105, = —1.08, n = —1 and for x > xq, where {xo,yo) is the 
elliptic fixed point. 

4 Dynamical properties of the Poinccire map 

Due to the the complexity of the perturbative term of the Hamiltonian of 
Eq.(7), the associated system in Eq.(5) is analytically intractable. In order 
to get some insight into its dynamical behaviour, we study the Poincare map 
(in fact the stroboscopic map) associated to this time-periodic Hamiltonian 
system. The unperturbed Hamiltonian HQ{r, 6) = r is globally degenerate, 
which means that d^Ho/dr^ = 0, Vr. The dynamical consequence of this 
degeneracy is that all the orbits are periodic having the same period. This is 
in contrast to KAM theory, where the unperturbed system is nondegenerate. 
Wc will show that after a periodic perturbation (periodic in time and in space) 
the corresponding Poincare map exhibits local and global bifurcations which in 
some cases are typical of the class of nontwist APMs which are perturbations 
of locally degenerate integrable APMs of an annulus, while in other cases are 
not. 



4-1 The nontwist property of the map 

Since the system given in Eq.(5) is time-periodic of period T = 2tt/uj, we 
define the function s : M — s(t) = t (mod 2n/u!), where is the circle 
identified with the interval [0, 27i/u). Hence, we have the autonomous system 
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of differential equations: 



X ^ y 

y — —X — e sin{ujx — us — Q;a;^/^sin us) (12) 
s = 1 . 



We denote by $j its flow. $j associates to each triplet (xq, |/o, Sq) G M x the 
position at the time moment t of the orbit starting at t — from [xq, yo, Sq). 
The plane 

E = {(a;,y,s) e X §i|s = 0} (13) 



is transversal to the flow and the map : E — >■ E deflned by P^{x,y) — 
$2,r/a)(-^' ^' 0) associated Poincare map. 

It is a classical result in dynamical systems theory that the Poincare map 
associated to a one degree of freedom periodic time-dependent Hamiltonian 
system is an APM. In order to interpret the dynamics displayed by the nu- 
merically computed Poincare map, let us deduce some geometrical properties 
of this map. The vector held F, associated to the system from Eq.(12), 

F[x, y, s) = {y, —x — esin(a;a; — lus — 6sin(a;s)), 1), (14) 



with (6 = aco^^^), is reversible with respect to the involution 7?. : x — > 
X S\ TZ{x, y, s) = {—x, y, 27r/a; — s), that is: 

FoTZ^ -DTZoF. (15) 



The reversibihty property (15) is equivalent to TZo^tt = $j o7^, Vt G R, which 
means that e R the diffeomorphism $j is 7?.-reversible, because = 

The Poincare map associated to the vector flcld F is i?-reversible, 
with respect to the involution i? : E ^ E, R{x,y) = {—x,y). Actually, notice 
that the reversor of the vector held F is related to the involution R by: 

n{x,y,s)^{R{x,y),27r/u-s). (16) 



It is clear that the deflnition of reversible maps is fulflUed by P^: 

Pe{R{x, y)) = y, 0) = ^2n/u>{T^{x, y, 27r/uj)) = ^^^^ 

n^-2^/^{x, y, 2t: /uj) = RP-'^{x, y). 
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Fig. 5. The nontwist property of the Poincare map. (a) The effect of the map on the 
semi-line defined by 6 = 7r/4; (b) The effect of the map on the semi-hne of ^ = 7r/2. 



The symmetry line useful in the analysis of the dynamical behaviour of the 
Poincare map is Fq = Fix(i?), having the equation a; = 0. 

The Poincare map Pq has the elliptic fixed point (0,0). After a slight per- 
turbation the fixed point persists as a symmetric elliptic fixed point of the 
reversible map Pg — it is a point (0, ^g) e Pq n Pi. In order to reveal the 
nontwist property of the map wc compute and visualize the effect of the 
map on different semi-lines emanating from the elliptic fixed point (Figure 5). 
As it results from Eq.(l), the rotation number for the map Pq is negative, 
more exactly p = —l/cu. Note that unlike the standard-like nontwist maps, 
here a line of constant 9 is not mapped onto a parabola, that is a curve with 
a single extremum, but onto an oscillating curve, having many minima and 
maxima. This fact will lead to a different behaviour of our nontwist map, in 
comparison to the dynamical properties of nontwist standard-like mappings. 

Besides this primary way to illustrate the nontwist property of the map Pg, we 
computed the self-rotation number of the orbits of the map P^, starting on the 
symmetry line Pq and using the turning angle method [16]. The self-rotation 
number for the map corresponding to the parameter e in the region of study 
is shown in Figure 6. Hence our map Pg is a nontwist APM, whose rotation 
number function p{y) is oscillating, having more than one extremum, that is, 
more than one invariant circle whose rotation number is a local extremum. 
Note that for the fixed parameter u = 3.0146 of the Hamiltonian system 
of Eq.(5), the associated Poincare map Pg has an elliptic fixed point which 
undergoes triplication. In Figure 7 we represent the triplication curve a*{e). 

In order to emphasize the significant differences in the phase portraits be- 
fore and after the triplication point, we add that in Figure 1 we illustrated 
the triplication bifurcation of the elliptic fixed point of the Poincare map for 
CO — 3.0146, e = 0.07 and a — —0.3 — panel (a) — a — —0.05 — panel (b) 
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Fig. 6. The self-rotation number of the orbits of points (0, y), with y G [—4,4]. 
The horizontal axis is the axis of y-values, while the vertical one corresponds to 
the self-rotation number. The dashed hne is at the level p = 1/3. (a) e = 0.01; (b) 
e = 0.02; (c) e = 0.03; (d) e = 0.04; (e) e = 0.05; (f) e = 0.06. 

— and a — 0.02 — panel (c). Prom the physical point of view it is only of 
interest the case in which a is constant and equal to ~ 0.3, and e varies. Thus, 
the considered pairs (e, a) are located above the triplication curve. According 
to [14] for this case a twistless circle can exist. Comparing the rotation num- 
ber associated to the cubic Henon map (Figures 3 and 4) and our Poincare 
map (Figure 6), we conclude that the Poincare map has an oscillating rotation 
number with a larger amplitude near the elliptic fixed point. This could be 
interpreted as a superposition of the nontwist property induced by triplica- 
tion on the nontwist behaviour due to the periodic perturbation of the initial 
Hamiltonian system. Near the elliptic fixed point, the Poincare map has a neg- 
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Fig. 7. The curve a*{e) whose points axe the triphcation thresholds of the central 
fixed point of the map P^. 

ative twist. The apparent discontinuities of the rotation number arc due to the 
distance chosen between the points on Fq whose rotation number is computed. 
While on the branch Tq = {{x, y)\x = 0,y > 0}, a pcriod-3 orbit has no point, 
it has one on the other branch, Fq = {{x,y)\x = 0,y < 0}. The orbit of a 
point on Fq encircling the homoclinic loops of the hyperbolic period-3 orbits 
can lead to an apparent discontinuity of the rotation number, although this 
function is continuous. The successive changes in the monotonicity of the map 
lead to a distinct nontwist behaviour in comparison with the known results 
reported for standard-like nontwist maps [4,9]. 

4-2 Birth of the dimerized island chains. Local and global bifurcations 

In our parametric study, we first considered the cases in which e is small, 
in order to study in detail the departure of the system from the simple har- 
monic oscillator as this parameter is increased. We restricted the study of the 
Poincare map to a rectangle limited by initial conditions close to physically 
reasonable values of the radius and velocity of the outer stellar layers. We 
focused on the creation of a typical phase-space structure of nontwist maps: a 
pair of dimerized island chains. This characteristic scenario of creation of new 
orbits in nontwist maps is illustrated in Figure 8. Panel a exhibits the phase 
portrait of the near-to-integrable map, where the elliptic fixed point is sur- 
rounded by invariant circles. At e ~ 0.03 the Poincare map has an invariant 
curve with three cusps (Figure 8b). Each cusp is a point of a new periodic 
orbit created through a saddle-center bifurcation for e ^ 0.04: a stable and 
an unstable fixed point emanate from each cusp (Figure 8c). The homoclinic 
loop of every hyperbolic point surrounds the corresponding elliptic one. Be- 
tween two successive hyperbolic points (in cyclic order) there is a heteroclinic 
connection. This separatrix structure called dimerized island chain appears in 
red in Figure 8c. In Figure 6d we have shown the self-rotation number p(0, y) 
for the case whose Poincare map appears in Figure 8c. In the former figure, 
the plateau p(0, y) = —1/3 corresponding to y E (1, 2) represents the constant 
rotation number due to the newly created elliptic fixed point on Fg . The 
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X X 

Fig. 8. The successive creation of the pair of dimerized island chains. Poincare maps 
for a = 0.3 and different values of e: {a) e = 0.01 — an "almost" harmonic oscillator; 
(b) e = 0.03 — the phase portrait just before the saddle-center bifurcation; (c) 
e = 0.04 — the creation of the first dimerized island chain; (d) e = 0.07 — the pair 
of period-3 dimerized island chain is complete. The homoclinic and heteroclinic arcs 
of the island chains are shown in red. 



dimerized island chain is also responsable for the minimum of p(0, y) around 

y ~ -2. 

Generically, dimerized island chains of the same period exist in pairs and are 
born in stages. The second period-3 dimerized island chain is created through 
the same bifurcation process at a higher value of e (Figure 8d). In Figure 6f, 
the small plateau of p{0,y) around |/ ~ — 1 explains the birth of the second 
dimerized island chain. Between the two island chains the invariant circles are 
meanders, that is, the radius along such a circle is not a univoque function of 
the angle. 

Around the last born period-3 dimerized island chain, a sequence of local 
and global bifurcations occurs. This is illustrated in Figure 9, where it can be 
seen that two independent orbits of the same period-3 are created by saddle- 
center bifurcation. They evolve in such a way that, finally, they interact with 
the orbits which belong to the first dimerized island chain. As e increases 
the newly born elliptic points approach the hyperbolic points of the dimerized 
chain. When e reaches a value of 0.11725 a global bifurcation occurs: the newly 
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Fig. 9. The global bifurcation: creation by saddle-center bifurcations of a pair of 
period-3 orbits for a = 0.3 and increasing values of e. The particular orbits involved 
in the bifurcations appear in red. (a) e = 0.0875 — the phase portrait at the 
threshold of a saddle-center bifurcation. It presents an invariant curve with six 
cusps; (b) e = 0.095 — two interwined orbits of period-3 have been created. Inside 
each new-born loop there is a period-3 elliptic point; (c) e = 0.11725 — the threshold 
of the global bifurcation. The topology of the separatrices has changed. 



created orbits interfere, and the hyperbolic points of the dimerized island chain 
become hyperbolic points with homoclinic eight-like orbits encircling the newly 
created elliptic points. 

In order to show the successive births of period-3 orbits as e increases with 
a = 0.3 fixed, we compute and visualize P^(Tq) and P^(Tq), k = 1,2,3. 
This process was first illustrated in Figure 5a where Tq and PKTq) = r2 
were represented. We show in Figure 10 the connection between the symme- 
try properties of the map and the creation of periodic orbits. A saddle-center 
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(a) (b) 




Fig. 10. The Poincare map before and after a saddle-center bifurcation. The fig- 
ures display also the symmetry lines F^^ (in blue) and F^^ (in green), with fc=0,l,2,3. 
(a) € = 0.055 — the phase portrait prior to the creation of the second dimerized 
island chain; (b) e = 0.07 — the second dimerized island chain has been created. 
The points satisfying {x, y) E F^ n F^ are period-3 points. 

bifurcation occurs at the value of the parameter e at which has a tangen- 
tial contact with F^. The corresponding Poincare map has for such an e an 
invariant curve with cusps which represent the points of tangency. Increasing 
the perturbation, the two symmetry lines intersect at two points, one being 
elliptic and the other, hyperbolic. 

4-3 Creation of chains of vortices 

A global bifurcation which is generic of nontwist maps is the so-called re- 
connection process, which is illustrated in Figure 11 for the nontwist multi- 
harmonic standard map [4]: 

Xn+i =Xn + ^-nuj + yl_^^ (mod 27r) (18) 
yn+i=yn + ksm[xn + arcsin(e siu x^)] , (19) 

with k = 0.2, e = 0.38 and different values of u. On the two sides of a twist- 
less circle, two Poincare-Birkhoff chains (necklaces of consecutive elliptic and 
hyperbolic points), having the same rotation number, approach each other 
(Figure 11a). At a given threshold, which is called the reconnection threshold, 
their hyperbolic points are connected by heteroclinic arcs (Figure lib). Vary- 
ing the parameter of the system further, two dimerized island chains emerge 
from the configuration created by the reconnection (Figure 11c). These chains 
are separated by meanders. 

In our case, we witness the creation of the nongeneric vortices or dipoles [9]. 
Figure 12 displays the birth of the first chain of vortices. As e increases, either 
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Fig. 11. Reconnection process in the nontwist standard- like multi-harmonic map. 
(a) 00 = 0.3287 — two independent Poincare-Birkhoff chains ; (b) u = 0.3297655 
— the reconnection threshold; (c) u> = 0.3315 — two dimcrizcd island chains. We 
illustrate also a meander separating the two dimerized island chains. 

Fg" or Fg intersects tangentially at a new point located at a radius larger than 
in the previous cases. As a consequence, new period-3 dimerized island chains 
are created after e crosses the value of tangential contact. However, each second 
dimerized island chain has a different creation mechanism. Unlike the previous 
case, the hyperbolic point of the first dimerized island chain bifurcates into two 
hyperbolic points in the transversal direction. In this process, a new elliptic 
point is born on the symmetry line. These hyperbolic points arc connected by 
three heteroclinic arcs: one surrounding the previously existing ehiptic point, 
another surrounding the new elliptic point, and the third one separating the 
two elliptic points. Thus, a pair of vortices is created. 

As e increases, the process of formation of chains of vortices continues for 
larger radii (Figure 13). Note that if the elliptic orbits of one chain of vortices 
intersect Fq", then the next one, which is created external to it, has a pair 
of elliptic points on the symmetry line Fq . The invariant rotational circle 
interpolating the hyperbolic points of the chain of vortices is the twistless 
circle. In Figure 13 it can be seen that such circles pass almost through the 
extremum of the symmetry lines Fg . As already said, the formation of vortices 
is not generic and therefore the twistless circle that appears between them is 
not either. 

4-4 Formation of the stochastic sea 

Our system is a clear example of weak chaos, where the perturbation itself 
creates a separatrix network at a certain eo (the dimerized island chains and 
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Fig. 12. Creation of the first chain of vortices for a = 0.3. (a) e = 0.0614 The 
saddle-center bifurcation has created the first dimerized island chain; (b) e = 0.06617 
— instead of a second dimerized island chain, a pair of vortices is born. See text for 
more details. 



the vortices chains) and then destroys it as e increases by producing regions 
of chaotic dynamics. The case of strong chaos implies that the unperturbed 
Hamiltonian intrinsically has separatrix structures and the perturbation just 
clothes them in thin stochastic layers. In both cases, the merging of the 
stochastic layers in the phase space may give rise to a single chaotic net- 
work called stochastic web. The larger e is, the wider is the stochastic web. 
Inside the cells of the web, there exists a set of islands of regular motion 
called invariant web-tori [17] whose dimensions are inversely proportional to 
the perturbation strength: for strong perturbations, they are engulfed by the 
stochastic sea. Our system presents this kind of behaviour, which is presented 
in Figure 14. 




Fig. 13. The phase portrait of the Poincare map corresponding to e = 0.2 and 
a = 0.3. One can see a sequence of concentric chains of vortices. The most exterior 
one is still incomplete. The figure displays also the symmetry lines F^, , /c=0,l,2,3. 
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Fig. 14. Formation of the stochastic sea(Q; = 0.3): (a) e = 0.3 — as the second 
dimerized island chain is destroyed, it is clothed in a stochastic layer; (b) e = 0.4 

— the stronger the perturbation, the wider the stochastic layers; (c) e = 0.5 — 
merging of the stochastic layers forming the stochastic sea. 



It has long been known (and not yet well understood) that the boundaries 
of the islands of regular motion can be "sticky" , meaning that chaotic orbits 
can spend a long time close to the islands and then escape in the chaotic sea. 
This phenomenon was observed already in [18] and was later called stickiness 
[19]. In general, it can be explained by the transformation of the web-tori into 
web-cantori as the perturbation is increased. Contrary to the absolute barriers 
constituted by the invariant tori, the cantori are only temporary barriers. That 
is, the orbits close to the main islands can reach the chaotic sea when they en- 
counter the gaps in the cantori, which happens after an indefinite time interval. 
Recent works [20,21] also suggest that the stickiness is due to the formation of 
higher-order resonant islands in the transition region between the main island 
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Fig. 15. A zoom of Figure 14b illustrating the islands-around-islands hierarchy. 



and the chaotic sea (Figure 15). The presence of these so-called dynamical 
traps appears to be generic for Hamiltonian systems and therefore their exis- 
tence critically determines the large-timescale behaviour of such systems. In 
this context, a behaviour extremely similar to sticky orbits has been shown to 
result from simulations of stellar variability using full hydrodynamical codes 
carried on significantly long time-intervals. In the absence of a better expla- 
nation, it was attributed to long-term (secular) nonlinear effects in the stellar 
envelopes. In a previous work [6], we have shown that our simple model recov- 
ers such a behaviour as it is generic of Hamiltonian systems. We also presented 
clear examples of sticky orbits together with an extensive discussion of their 
implications in the framework of the classification of variable stars. Therefore, 
we mention here only the fact that if a Hamiltonian approach is intrinsic to 
the phenomenon of stellar variability, we can expect that different categories 
of variable stars (classified according to their amplitude, frequency, irregular- 
ity) may be viewed under a more unifying perspective than it is presently 
considered. Moreover, considering that the sporadic excursions of the sticky 
orbits to the chaotic sea translate into high velocities of the outer layers, mass 
loss is very likely to occur, again in accordance with the observations. 

Another important and distinctive feature of the phase portrait of nontwist 
maps is the existence of meanders, as it was mentioned before when discussing 
Figure 8d. Meanders are created between two successively born dimerized 
island chains or between two chains of vortices. In nontwist standard-like maps 
they become usual invariant curves after the reconnection of the two chains. 
For the sake of clarity, we illustrate next the relationship between the meanders 
and the reconnection process for our map. In Figure 16 we show the creation of 
meanders from the reconnection of the Poincarc-BirkhofT chains of period 34. 
In the case of the standard nontwist map [9], in the reconnection of Poincare- 
Birkhoff chains of even periodic orbits, periodic points of the chains having 
the same stability type are aligned in phase. In other words, to an elliptic 
(hyperbolic) point of one chain corresponds below or above also an elliptic 
(hyperbolic) point. Due to this alignment, the reconnection is nongcncric and 
corresponds to a hyperbolic-hyperbolic collision leading to the formation of 
vortices. In our case, the two Poincare-Birkhoff chains of even period (period 
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Fig. 16. The reconnection process: (a) e = 0.405 - The Poincare-Birkhoff chains 
confirming the generic process of reconnection illustrated in Figure 11; (b) e = 0.4 
— After the reconnection, the resultant meander with 34 cusps, slightly visible in 
Figure 14b. 



34) have periodic points of opposite stability type which approach each other in 
their way to reconnection (Figure 16a). Therefore, the reconnection is generic 
and in the subsequent dimerized islands chains, hyperbolic-elliptic collision 
occurs. Moreover, after the reconnection, the twistless circle turns from a graph 
of a function of the angular variable to a meander (Figure 16b). This meander 
is slightly visible in Figure 14 b. 

Meanders appear to be robust invariant circles acting also as separatrices be- 
tween distinct regions of chaotic dynamics. For stronger perturbations, their 
destruction allows the chaotic orbits to reach other stochastic regions pre- 
viously separated by the meander. This behaviour was observed in nontwist 
standard-like maps [10], but until now there is no explanation for this robust- 
ness. 



5 Conclusions 



In the present work we have extended the preliminary results of Ref. [6] con- 
cerning the dynamics of a forced oscillator as a model of irregular stellar 
pulsations. The driving is characterized by two parameters, the fractional am- 
plitude of the internal perturbation a and the total amplitude of the driving 
e. Our aim was to undertake a more detailed analysis of the transition to the 
nontwist property in the associated Poincare map. We have demonstrated that 
for a given e > there exists an a*{e) that corresponds to a triplication of the 
central elliptic fixed point of the map. As we were interested in values of the 
parameters a and e above the triplication curve, the nontwist character of the 
map was undoubtly present [14] . 

Special attention was devoted to the process of formation of periodic orbits 
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and to periodic-orbits collision. The reversibility of the Poincarc map and 
consequently, the symmetry properties allowed a clear identification of the 
structure of the periodic points. More precisely, the evolution and the inter- 
sections of the symmetry fines of the system provided good tracers of the steps 
of creation of dimerized island chains and vortex chains. We strongly argued in 
favor of these particularities of the map being entirely due to the superposition 
of two factors: the nonmonotonicity of the rotation number function induced 
by the triplication of the elliptic fixed point and the nonmonotonicity of the 
same function due to the oscillating perturbation. The latter feature admits 
alternative occurrence of maxima and minima, where the twist condition is 
violated. This translates into the minima and maxima of the symmetry lines 
which transform the invariant curves passing through them into twistless cir- 
cles. Using these elements, we have followed the local and global bifurcations 
until the formation of stochastic layers around separatrices together with the 
associated sticky orbits. We have proved that the reconnection of even peri- 
odic orbits is generic by following the reconnection of two Poincare-Birkhoff 
chains of period 34. From the astrophysical point of view, we have dwelled on 
the implications of the sticky regions and of the interplay between stochastic 
and regular regions in the framework of the analysis of stellar variability. In 
order to provide quantitative characterization of sticky orbits in the frame- 
work of stellar variability, we shall concentrate in the future on the causes and 
probabilistic estimates of the stickiness mechanism. 
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